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Abstract 

We report the analysis of the first superburst from a transiently accreting neutron star system with the All-Sky Monitor (ASM) 
on the Rossi X-ray Timing Explorer. The superburst occurred 55 days after the onset of an accretion outburst in 4U 1608-522. 
During that time interval, the accretion rate was at least 7% of the Eddington limit. The peak flux of the superburst is 22 to 45% 
of the Eddington limit, and its radiation energy output is between 4 ■ lO'*' and 9 • 10^' erg for a distance of 3.2 kpc. Fits of cooling 
models to the superburst light curve indicate an ignition column depth between 1.5 • 10 and 4.1 • lO'^ g cm^^. Extrapolating the 
accretion history observed by the ASM, we derive that this column was accreted over a period of 26 to 72 years. The superburst 
characteristics are consistent with those seen in other superbursting low-mass X-ray binaries. However, the transient nature of the 
hosting binary presents significant challenges for superburst theory, requiring additional ingredients for the models. The carbon that 
fuels the superburst is thought to be produced mostly during the accretion outbursts and destroyed in the frequent type-I X-ray bursts. 
Mixing and sedimentation of the elements in the neutron star envelope may significantly influence the balance between the creation 
and destruction of carbon. Furthermore, predictions for the temperature of the neutron star crust fail to reach the values required for 
the ignition of carbon at the inferred column depth. 
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views, see lLewin et alll 19931 and IStrohmayer & Bildsten"2006l). 
The flash layers heat up during a fraction of a second and cool 
down in an exponential-like manner lasting 10 s to a few min- 
utes. As a result, a spectrum may be observed which is typ- 
ically well fit by a black body with a temperature decreasing 
from a peak of typically kr = 2 to 3 keV and an emission re- 
gion similar in size to what is expected for a NS dSwank et al.l 
19771) . Some nuclear energy generation may persist throughout 
the cooling phase, prolonging the bursts somewhat. A number 
of bursts exhibit photospheric radius expansion (PRE) where 
the burst flux is thought to reach the Eddington limit. From the 
observ ation of PRE bursts from 4U 1608-522. iGallowav et al.l 
(l2006h derived a distance to this source of 3.2 ± 0.3 kpc. During 
5 PRE bursts, burst oscillations were detected at 619 Hz, which 
implies the fastest known spin frequency for accreting neu- 
tron stars (Gallowav et al. 2006), with the possible exception of 
XTE J1739-285 tKaaret et ah,2007i) . 



Since a few years a different kind of type-1 X-ray burst is 
detected in abo u t 10% o f all bursters: so-ca lled 'superbursts' 
dCornehsse et alj|2000t IStrohmaver & Brownil2002D . These are 
about 10^ times more energetic and last for hours to a day. 
The y are thought to restilt from flashes of carbon-rich lay - 
ers ( IStrohmaver & Brownl 120021: ICumming & BfldstenI l200lk 
These layers are presumed to be much further down in the 
neutron star than those where ordinary X-ray bursts occur, 
which would explain the longer duration. It has been ob- 



1. Introduction 

4U 1608-522 is a bright lo w-mass X-ray binary (LMXB) with a 



neutron-star (NS) primary ( Behan et al. 1976: Tananbaum et aXl 
19761) . QX Nor was identified bv lGrindlav & LiUeii (1 1 9781) as flie 



optica l counterpart. The acc retion onto the neutron star is tran- 
sient (iLochner & Roussel-D upre 1994 and references therein) 
and the outburst behavior complex. Three mai n flux states can be 
resolved in 4U 1608-522 dWachter et al.l2002l) : an outburst state, 
a low-intensity state and a true quiescent state. In contrast to 
many other LMXB transients, its true quiescent state has a rela- 
tively low duty cycle of only about 50%. Wachter et al. speculate 
that the mass transfer rate might be fluctuating around the critical 
threshold separating stable accretion systems from true transient 
systems, possibly due to varying stellar spots on the donor star or 
to a varying vertical accretion disk structure which may change 
the hi'adiation and, th us, the mass transfer rate through the disk. 
IWachter et al.l (|2002|) also identified a modulation in the optical 
light curve with a period of 12.9 hr which they suspect to be near 
the orbital period. Thus, the donor is probably hydrogen-rich. 

4U 1608-522 exhibits so-calle d type-1 X-ray bursts 

dBehan et al.l 119761 : iTananbaum et al.l Il976i) . These are under- 
stood to result from thermonuclear shell flashes on the sur- 
faces of NSs of layer s of accreted hy drogen and helium 
dWooslev & TaamI 119761: iMaraschi feCa vahere 1977; for re- 
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served that superbursts influence the normal bursting behavior. 
Each time when the start of a superburst has been observed 
with sufficient statist ics a precursor burst was seen (see e.g. 
IStrohmaver & Brownl 12002). Furthermore, after the superburst 
the normal burs ting behavio r is quenched for approximately one 
month (see e.g. lKuuUcers et al. 2002b). 

Superbursts have so far only been seen in systems in which 
the NS has been continuo usly accreting for at least 10 yr 
jKuulkersI 120041; lin 'tZand e t al. 2004), although not in every 
such system ( Keek et alT 2006i) . Here we present a superburst 
from 4U 1608-522, which represents the first case for a 'clas- 
sical' transient. A superburst has also been observe d from the 
transient source KS 1731-260 dKuulkers et al.l2002bl) . However, 
unUke KS 1731-260, 4U 1608-522 exhibits states of high flux 
with a duration that is short with respect to the expected super- 
burst recurrence time. This may have interesting consequences 
for the theory of superbursts, because the different states of the 
transient allow us to study different levels of carbon production 
and burning in one source. Perhaps most importantly, whereas 
the crust in KS 1731-260 is expected to have been heated signif- 
icantly out of thermal equilibrium with the core during its long 
outburst, the crust of 4U 1608-522 is not expected to be signif- 
icantly heated during its short outbursts. This presents a signifi- 
cant challenge to carbon ignition models for superbursts, which 
require a hot crust to achieve ignition at the depths inferred from 
observations of superbursts. Despite these differences, the su- 
perburst from 4U 1608-522 has similar properties to superbursts 
from other sources. 

The paper is organized as follows. First we describe the ob- 
servations and spectral calibration in Sect.|2] In Sect.|3]we report 
on an analysis of the X-ray burst behavior of 4U 1608-522 as a 
function of its apparent mass accretion rate. Next, in Sect. |4] we 
present the analysis of the detected superburst. Finally, we dis- 
cuss the implications of the results for superburs t theory in Sect. 
[5l Init ial r eports of this sup erburst appeared in iRemillard et al.l 
(l2005h and lKuulkersl(l2005l) . 



2. Observations and spectral calibration 

2.1. Observations 

To study both the long-term accretion and bursting behavior of 
4U 1608-522, as well as the superburst, we employ observations 
performed with multiple instruments on-board four X-ray obser- 
vatories. 

The All-Sky Monitor (ASM) on the Rossi X-ray Timing 
Explorer (RXTE) jLevine et al J 19961) consists of three Scanning 
Shadow Cameras (SSC), each containing a position-sensitive 
proportional counter The cameras are mounted on a rotating 
drive. Data are accumulated in so-called dwells of 90 seconds. 
After each dwell the rotation drive changes the orientation of 
the SSCs. Spectral information is available from three channels 
with corresponding energy ranges of roughly 1.5-3, 3-5 and 
5-12 keV. After the launch of RXTE in December 1996 until 
November 2006, the ASM observed 4U 1608-522 for a total of 
2.8 Ms. On May 5th 2005 the ASM observed the superburst from 
4U 1608-522. 

RXTE also carries the Proportional Counter Array (PCA). 
The PCA consists of five proportional counter units (PCUs) with 
a total geometric area of 8000 cm^ and has a bandpass of 1 
to 60 keV CJahoda et al. 2006 ). 4U 1608-522 was observed for 
1 .6 Ms in total, mostly when it was out of qu iescence (e.g., 
Ivan Straaten et al ] |2003U Gierlinski & Donell2002h . 



The BeppoSAX observatory (iBoeUa et al.l Il997ah was 
launched in April 1996 and carried two Wide Field Camera's 
(WFCs) as well as four N arrow Field Instruments (NFI). The 
WFCs dJager et al.l 1 19971) are coded mask aperture cameras 
with a band pass of 2-28 keV. During the six year lifespan 
of BeppoSAX a campaign of semi-yearly observations of the 
Galactic Center were carried out, which resulted in an exposure 
time of 3.8 Ms for4U 1608-522. 

The BeppoSAX NFI performed two pointed observations 
of 4U 1608-522. The first observation was performed dur- 
ing an outburst on February 28 1998. In this paper we 
analyze the broad-band spectrum obtained with the follow- 
ing NFI in the indicated energy bands: the Low Energy 
Concentrator Spectrometer (LECS; jParmar et all 119971: 0.12- 
4 ke V), the Mediurn Energy Concentrator Spectrometer (MFCS; 
Boe lla etal.lll997bl: 1.8-10 keV) and the Phoswich Detection 
System (PDS: lFrontera et al.|[T997 : 15-220 keV). The exposure 
time depends on the instrument and is 30 ks for the MFCS. The 
second observation took place when 4U 1608-522 was in quies- 
cence. Since we are interested in the production of carbon from 
accreted matter, we forgo the analysis of the data from this ob- 
servation. 

On May 5th 2005, 33 minutes before the start of the 
ASM observation of t he superburst, the High Energy Transient 
Explorer 2 (HETE-2: IRicker et alJl2003l) observed a flare from 
4U 1608-522 with two instruments : the Wide field X-ray 
Monitor (WXM: IShirasaki et alj |2003) with a 2 to 25 keV band 
pass a nd the French Gamma Telescope (FREGATE: lAtteia et al.l 
l 2003l) which has a band pass of 6 to 400 keV. 

The Inter national Gamma-Ray Astrophysics Laboratory 
(INTEGRAL; IWinkler et al. ,2003.) was launched in October 
2002. The Imager on Board the Integral Satellite (IBIS; 
lUbertini et al.1 120031) is a coded aperture camera containing 
two detectors. We use results obta ined with the INTEGRAL 
Soft Gamma-Ray Imager (ISGRI; iLebrun et all 120031) detec- 
tor, which has an energy range from 15 to about 500 keV. 
IBIS/ISGRI has observed 4U 1608-522 up to September 2005 
for a total of 6.3 Ms. 

In this paper we analy ze spectra usin g version 1 1 .2.0 of the 
XSPEC software package dArnaudl 19961) . 

2.2. Spectral calibration of the ASM 

The spectral response of the ASM detectors is not well defined. 
A pre-flight effective area array is av ailable, while a redistribu- 
tion matrix is not (see discussion in 'Kuulkers' '2002") . We con- 
struct such a matrix by modeling, for an infalling photon with 
a certain energy, the energy distribution over the three channels 
by a Gaussian with a full width at half maximum of 20% of the 
centroid energy, which is a typical wi dth for the type of propor- 
tional counter used in the ASM (e.g. Fraser"l989^. The response 
matrix thus obtained is merely a first order estimate, but suffices 
for rough calculations. 

We employ the Crab source to investigate the accuracy of 
our response matrix. The average count rate over all ASM ob- 
servations from this source in the full 1.5-12 keV band-pass is 
75 .4 c s^ ' . The X-ray spectrum can be described by an absorbed 
power law. iKirsch et al.l (l2005b performed simultaneous model 
fits to Crab spectra obtained with 22 X-ray instruments. Using 
the results of these fits in the 2-10 keV range, our response 
matri?|3 predicts an ASM 1.5-12 keV count rate of 70.5cs^'. 



^ iKuuIkersI ( |2002|) finds the count rate is over-predicted. This is due 
to a diiJerent choice of values for the spectral model parameters. 
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Figure 1. 1.5-12 keV RXTE-ASM light curve of 4U 1608- 
522 at a 2-week time resolution. Data points with errors in ex- 
cess of 0.5 c s^' were excluded from this plot. The connect- 
ing line is broken if data points are more than 2 weeks apart. 
The vertical lines indicate 37 bursts detected with the WFC, 19 
with the ASM (dashed lines ind icate tentative bursts), 31 with 
the PCA (iGallowav et all |2006|) and 19 with the IBIS/ISGRI 
jChelovekov et al. 2007 ). The long vertical line indicates the 
time of the superburst. 'BSAX' indicates the times of two 
BeppoSAX observations. 



3. History of accretion and X-ray burst activity 

3.1. Long-term light curve 

Figurefflshows the 1 .5-12 keV ASM Ught curve of 4U 1608-522 
Clearly visible are the three states identified by IWachter et alj 
(2002; see Sect. [Til. Four major outbursts, where the peak photon 
count rates was in excess of 20 c s^' for more than one month, 
are visible in the 1 1 -year time span of the observations (see Sect. 
l3.2l for a discussion on the selection criterion). The first outburst 
was ongoing at the start of the ASM observations and lasted 
until 60 days afterward, with fluxes in excess of 20c s^'. The 
following outbursts lasted with fluxes above 20 c s^' for 44 d 
(MJD 50848-50892), 48 d (MJD 52475-52523), and 80 d (MJD 
53438-53518), respectively. There are also five minor outbursts, 
after the two latter major outbursts. The low-intensity states are 
visible in the first half of the mission. From comparing the first 
half and the latter half of the light curve, it appears that either 
a low-intensity state emerges after a major outburst, or a series 
of minor outbursts. This characteristic is only now apparent, af- 
ter 1 1 years of observation (Wachter et al. 2002 only considered 
the first 5 years of the data set). We determined the average flux 
over the complete ASM data to be 2.466 ± 0.006cs"'. This is 
about 2% of the flux reached during the brightest bursts seen in 
the ASM. 



3.2. Accretion rate 



Therefore, when performing spectral analyses using this matrix, 
the normalization of the models needs to be corrected by a fac- 
tor 1 .07 to account for this discrepancy. Note that in principle 
this factor can be different for each of the three SSCs and can 
vary with time. However, we find that in a time interval of 100 
days centered at the start of the superburst the difference in Crab 
count rate between any two SSCs never exceeds the 3a level. 

Apart from the accuracy in predicting the count rate of the 
Crab, we also investigate how well our response matrix can 
reproduce the typical model parameters that are found by fit- 
ting an absorbed power law to the three-channel spectral data. 
We extracted a spectrum from all the ASM data on the Crab 
available at the time of writing. We fix the hydrogen absorp- 
tion column den sity at t he val ue of A^h ~ 0.45 • 10^^ cm^^, 
as found by Kirsch et all (l2005h (in the 0.1-1000 keV energy 
range). FollowinglKirsch et all we use the abundances found b y 
IWilms etalJ (12000) ' and cross sections from lVerner eTal] (Il996l) . 
Leaving free the power law index and the normalization, we do 
not find an acceptable agreement with the data. Only if we add in 
quadrature 10% of the flux to the uncertainty of each data point 
are we able to obtain an acceptable fit with x^^d — ^ ■ Taking into 
account the correction factor derived previously, we find the best 
fit with a photon index of F = 2.01 ± 0.12 and a normalization 
of A^poweriaw = 9 ± 2 photo ns keV~^cm~^s~^ which is consis- 
tent with the results from Kirsch et al.l (l2005h . The uncertainty 
in A'poweriaw is large, because by definition A^poweriaw is flie pho- 
ton flux at 1 keV, which is outside of the ASM energy range. The 
XSPEC power law model 'pegpwrlw' uses a user-defined energy 
range for the normalization. Employing the ASM bandpass gives 
an uncertainty in the normalization of 6%. 



When matter is accreted onto the neutron star surface, some of 
the accretion energy and/or some of the material may leave in a 
jet. If these losses are significant, the persistent flux is not a good 
trace r of the mass accretion rate M. For black-hole X-ray bina- 
ries, iFender et all (l2005h argue that jets are not present during 
the high soft state. It is during the high soft state of 4U 1608- 
522 that we are particularly interested in the precise value of the 
accretion rate. Therefore, assuming that the same holds true for 
neutron star binaries, in this paper we will assume there are no 
losses through a jet and that the persistent flux is a good measure 
for M. Expressed as a fraction of the Eddington-limited mass ac- 
cretion rate Medd- which is MEdd = 2 • lO^'^Mgyr^' for a canon- 
ical hydrogen-accreting neutron star with a mass of 1.4M0, it 
can simply be derived as M/MEdd — F /pEddi with F the bolo- 
metric flux and fedd the Eddington-limited flux. We will deter- 
mine the bolometric flux below. The Eddington-limited flux is 
exhibited during photospheric radius e xpansion (PRE) bursts . 
RXTE PCA observed 12 such bursts dGallowav et all l2006l) . 
The mean unabsorbed bolometric peak flux of these bursts is 
FEdd= (1.32±0.14)10-^ergs-icm-2. 

Van Straaten et al. (l2003l) make an extensive study of the 
spectral changes of 4U 1608-522 during the various states and 
the associated timing properties. They conclude that 4U 1608- 
522 is an Atoll source that exhibits the 'banana' state during out- 
bursts and moves in the color-color diagram to the 'extreme is- 
land' state when the flux goes to levels of the low-intensity state. 
The distinguishing parameter is the 'hard color' (i.e., the flux be- 
tween 9.7 and 16.0 keV divided by the flux between 6.0 and 9.7 
keV) which is about 1.15 in the extreme island state and 0.55 
in the banana state. This implies that one-band measurements, 
such as the ASM count rates, are not an accurate measure of the 
bolometric flux, let alone the accretion rate. In order to derive 
the bolometric flux during the low and high accretion states, we 
analyze the broad-band outburst-observation with BeppoSAX as 
well as a large number of observations with the RXTE PCA. 
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Table 1. Parameters and 1<7 uncertainties of fitting the out- 
burst spectrum of 4U 1608-522 from BeppoSAX NFI observa- 
tions on 28 Feb 1998 with an absorbed disk black body" and 
comptonized'' model. 



A'bb 
kTo 

X 

Atomptt 

Flux (0.01-300 keVf 
ASM count rate'' 



(8.91 ±0.05) 10^' cm" 

2.38 ±0.07 keV 
18.3±0.9 

0.478 ±0.005 keV 
3.6±0.2keV 
3.7±0.2 
0.40 ±0.02 



1.40/160 
(1.75±0.10) 10^ 
39.4±0.9cs-i 



' erg s ' cm ^ 



" Disk black body model parameters: krBB is the temperature at the 
inner disk radius R^^ and A'bb = ('^ki/'^fo) ^ normalization, with 
Rj„ in units of km, d\Q the source distance in units of 10 kpc and 6 the 
inclination. 

Comptonized model parameters: kTo, kr^ are resp. the seed photon 
and plasma temperature, T the plasma optical depth for a disk geometry 
and A'comptt the normalization. 

Unabsorbed 0.01-300 keV flux, which we take to be bolometric. 

Average ASM count rate during the observation. 



Two broad-band observations were carried out during the 
ASM coverage with BeppoSAX, one of which during outburst 
(for the timing of these observations, see Fig. [U. To deter- 
mine the flux during this observation we analyze spectra ob- 
tained with the LEGS, MFCS and PDS. We subtract background 
spectra taken from blank field observations at the same detec- 
tor positions for the LFCS and MFCS spectra, while for the 
PDS off-source pointings are used as background. The spectra 
are rebinned to obtain at least 15 photons in each bin while 
sampling the instrument resolution with at most 3 channels 
at all energies, to ensure the applicability of the statistic. 
An error of 1% is added in quadrature to the statistical error 
per bin to account for systematic uncert ainties. This is com- 
mon practice in BeppoSAX analyses (Fiore et aP Il999h . We 
fit the data from all three instru ments simultaneou sly with the 
generic LMXB model (see e.g. ISidoh et akl 20011). consisting 
of a multi-temperature disk black body (iMitsuda et al.l Il984t 
iMakishima etal.1 119861; 'disk bb' in XSPEC) in combination 
with a comptonized spectrum (Titarchuk 1994; Hua & Ti tarchukI 
Il995t iTitarchuk & Lyubarskij 1995; 'comptt' in XSPFC), both 
absorbed by cold interstellar matter fol lowing the model by 
iBalucinska-Church & McCammonl (|1992|) Cphabs' in XSPFC). 
During the spectral analysis we allow the normalization between 
the three instruments to vary and find that the best fit values are 
within acceptable limits; the normalization between the LFCS 
and MFCS is 0.897 ±0.003, while 0.7 to 1.0 is acceptable; 
between the PDS and MFCS the normalization is 1.21 ±0.04, 
while 1.1 to 1.3 is acceptable (Fiore et al. 1999). The results of 
the fit as well as the derived flux are provided in Table [T] The 
flux is determined by extrapolating the model from the observed 
range to 0.01-300 keV. Extrapolating the model further while 
correcting for the interstellar absorption does not lead to a sig- 
nificant increase in flux. Under the assumption that no other fea- 
tures in the spectrum other than those we modeled have a sub- 
stantial flux contribution, we refer to this as the bolometric flux. 



From the average ASM flux during this observation, we find that 
1 ASM c s"' is equivalent to (4.4±0.2) • lO^^'erg s^'cm"^. 

To further study the bolometric correction that may be ap- 
plied to the ASM data as a function of the persistent flux, we 
employ 160 observations carried out with the RXTE PCA be- 
tween March 1996 and February 2006. The PCA spectra are 
extracted from the standard products, with exposure times be- 
tween 96 s and 13 ks and a total exposure of 310 ks. They are 
modeled between 3 and 40 keV (the calibrated bandpass) with a 
generic model consisting of a black body a nd a cutoff power law, 
both absorbed following the model by Mo rrison & McCammonl 
([19831) with A^H = 8 • 10^' cm"^. The resulting fits are accept- 
able in 144 cases with x^^d < 2; in the remaining cases x^^i rises 
to up to 6. Despite the sometimes bad fit, we employ all ob- 
servations since we are interested only in obtaining constraints 
for the bolometric correction and not in physical implications of 
the spectral model. The fit outcome is tested leaving A^h free; 
no evidence is found for fluctuations in A^h above 5 • 10^^ cm^^. 
We determine the photon flux between 1.5 and 12 keV before 
correcting for absorption and normalize it to observed ASM 
fluxes through calibration with Crab fluxes and by cross corre- 
lating the PCA measurements for 4U 1608-522 with contempo- 
raneous ASM fluxes. We calculate the 0.1 to 100 keV energy 
flux after a correction for absorption. Obviously, this implies a 
fairly large extrapolation from the 3 to 40 keV PCA bandpass 
and we check how much flux is found outside that bandpass. It 
varies between 10% at the brightest state to 60% in the faintest 
state. This indicates what the maximum error is of the extrap- 
olation. The results are plotted in Fig. [2l I t clea rly shows the 
two states identified by Ivan Straaten et al.l (l2003b . The bound- 
ary between both states is at 5 to 10 ASM c s^^ At the bright 
state, when the spectrum is relatively soft, I ASM c s^Ms equiv- 
alent to 5.5 • 10^"^ erg s^'cm^^ and at the hard faint state to 
1.0-10^'^ erg s^'cm^^. For the bright state we can calibrate this 
result using the BeppoSAX broadband observation of the 1998 
outburst from which we derived above the equivalence of 1 ASM 
c s^' to (4.4 ±0.2) • lO^"' erg s^'cm^^. We use this number as 
conversion factor for the high flux state. 

Applying the bolometric correction factors to the ASM 
flux, we find that the average persistent flux was (2.313 ± 
0.013)10^^ergcm^^s^\ while the peak flux reached during 
the three large outbursts was 2.2 ■ 10^** to 3.0 • 10^** erg cm^^ 
s^'. The mass accretion rate is on average M/Mem = (1-8 ± 
0.2)10"^ and 0.17 to 0.23 in the peaks of the outbursts. 

Using the conversion factor of ASM counts to bolometric 
flux for count rates above 10 c s^' we determine the fluence of 
the three major outbursts observed by the ASM for the duration 
that the count rate exceeded 20 c s^'. For the last outburst we 
also determined the fluence up to the moment that the superburst 
occurred. 20 c s^^ corresponds to an accretion rate of M/Mgdci = 
(6.7 ±0.8) 10^2 and is chosen as a lower limit since it is close 
to the value of 0. 1 that is required for the production of carbon 
(ICumming & Bildsten 2001), while selecting the periods of high 
flux in one interval for each outburst. The bolometric fluences 
as well as the peak fluxes during the outbursts are presented in 
Tabled 



3.3. Long-term bursting behavior 

Indicated in Fig. [1] are the bursts detected from 4U 1608-522 
with the WFC on BeppoSAX until the mission ended in April 
2002, the IBIS/ISGRI o n INTEGR AL until September 2005 
dChelovekov et alJl2007l) . the PCA dGalloway et alJl2006l) and 
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Outburst 


Start time^ 


Duration^ 


Fluence 


Energy'' 


Peak Flux 


PeakM/MEdd'' 




(MJD) 


(days) 


(lO^^ergcm"^) 


(lO^^erg) 


(10"**ergs"'cm"2) 




1 


50848.0 


44.0 


6.3(3) 


0.77(15) 


3.0(2) 


0.23(3) 


2 


52475.0 


48.0 


6.9(3) 


0.85(16) 


2.6(2) 


0.20(3) 


3 


53437.5 


79.5 


10.4(5) 


1.3(2) 


2.24(19) 


0.17(2) 


3 before SB 


53437.5 


57.6 


8.2(4) 


1.01(19) 


2.24(19) 


0.17(2) 



^ Start time and duration of time interval when the outburst flux exceeds 20 ASM c s^ . 
Energy calculated from fluence using distance d = 3.2 ±0.3 kpc. 

Unabsorbed bolometric peak flux determined from the 1-day ASM light curve (10"'*ergs"'cm"2). 
Peak mass accretion rate in units of the Eddington limit assuming a 1 .4 Mq neutron star. 




1 10 

ASM flux [c s"'] 

Figure 2. Flux as a function of the ASM count rate. The distinc- 
tion between the two spectral states is clearly visible. The two 
lines are not fits, but act as guides to distinguish the two states. 



the ASM on RXTE. One burst was detected by both the WFC 
and PCA. The ASM bursts are identified through 7 ) searching 
individual 90-sec average dwell data of 4U 1608-522 for iso- 
lated high points that are at least 3(7 above a 100-day running 
average, resulting in roughly 200 burst candidates; 2) studying 
the raw 1-sec resolution light curves of the candidate dwells, to 
eliminate candidates that have no clear flare feature with a time 
scale smaller than the 90-s dwell time; 3) searching for flares 
with clear fast-rise exponential-decay profiles, resulting in the 
identification of 7 certain bursts. The remaining 12 candidates 
are classified as tentative bursts from 4U 1608-522. 

For studying the bursting behavior as a function of the per- 
sistent flux, Fig.Oshows a histogram of the ASM-measured per- 
sistent flux when a burst is seen for all bursts observed with the 
WFC, PCA and IBIS/ISGRI as well as the certain bursts ob- 
served with the ASM. There is the suggestion that bursts are 
rare during major outbursts while much more frequent during 
the extended low-intensity states during MJD 50150-50500 and 
50897-51800. Quantitatively, however, the data are inconclusive 
in this respect. The PCA data set is somewhat biased, because 
those observations were targeted at 4U 1608-522 mainly when 
it was in outburst. The other data sets are serendipitous and, 
therefore, present an objective view of the bursting behavior The 
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Figure 3. Number of observed bursts (solid line) as a function of 
persistent flux Fpeis ■ For the latter we use the average ASM count 
rate in a one-day interval around the moment the burst occurred 
converted to bolometric flux. This includes all bursts observed 
with BeppoSAX WFC, RXTE PCA, INTEGRAL IBIS/ISGRI 
as well as the clear bursts identified in the RXTE ASM light 
curve. The dashed lin es indicate the region of Fpers where 
iRevnivtsev et^ (120011) observed mHz QPOs (see Sect. ISH l. 
The dotted line represents the total PCA, WFC, IBIS/ISGRI and 
ASM exposure time as a function of the ASM count rate. For 
PCA, WFC and IBIS/ISGRI observations we use the average 
ASM count rate during an observation. Together with the burst 
histogram, one observes that the burst rate does not vary signif- 
icantly as a function of accretion rate in the range where type-I 
bursts occur 



WFC exposure time during major outbursts is 21.3 times smaller 
than during low-intensity states, while the ratio in the number of 
bursts is 1 to 36. 2 out of 19 bursts observed by IBIS/ISGRI take 
place during a major outburst, while the source is observed in 
outburst during 12.7% of the exposure time. The ASM data are 
also inconclusive: 2 out of 19 bursts occurred at 1-day average 
fluxes of 10 c s^' while the fraction of the total time spent in 
that flux regime is 7.4%. The WFC, PCA, IBIS/ISGRI and ASM 
light curves exhibit bursts with an average recurrence time of, 
respectively, 1.19, 0.60, 3.79 and 1.71 days. 

The bursting behavior is characterized by a: the ratio of 
persistent fluence between bursts and the burst fluence. Using 
the PCA observations we calculate the average value of a both 
for the outbursts, when the ASM count rate was in excess of 
20 cs^^ and for all observations when the ASM count rate 
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Figure 4. Superburst light curves where each data point repre- 
sents a 90s dwell with its la uncertainty, a) 1.5-12 keV ASM 
count rate at a 0.5 day time resolution. Shown is the outburst 
in which the superburst takes place, b) Zoomed-in part of a. 
Indicated are intervals I and II of the persistent flux prior to the 
burst, c) hardness ratio for each dwell, defined as the ratio of the 
counts in the 5-12 keV and the 1.5-3 keV energy bands. One 
data point with an error in excess of 4 is excluded from the plot. 
d) Zoomed-in part of top figure on the superburst. Burst intervals 
1-4 are indicated as well as the exponential decay fitted to the 
first 5 hours of the burst with an e-folding decay time of 4.8 hr 
(see TableO. 



was be low this level, but not during quiescence. iGallowav et alj 
('2006) determined the fluence of 3 1 of the type-I X-ray bursts 
observed with the PCA. We find the persistent fluence between 
bursts by multiplying the average ASM flux during the PCA ob- 
servations by the burst recurrence time. For the latter we take the 
ratio of the total PCA exposure time and the number of observed 
bursts during either the outbursts or the low-flux state. During 
the outbursts the PCA observed only one burst. Therefore, from 
the PCA data we only have a lower limit to the burst recur- 
rence time in the outbursts, from which we find the lower limit 
a > 745 ± 9. The PCA observations when the ASM flux was 
less than 20 c s^' yield an average value of a = 51.5 ±0.2. 

4. Superburst analysis 

4.1. Light curve 

In Fig. m we show the ASM light curve of 4U 1608-522 for 
the outburst in which the superburst takes place. The onset of 
the superburst is not observed by the ASM, as it falls in a 1.5 
hour data gap. Thirty-three minutes before the start of the su- 
perburst observation by the ASM, the light curves of the WXM 
and FREGATE onboard HETE-2 exhibit a fast rise of 3.6 s fol- 
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Figure 5. HETE-2 light curve of likely onset superburst. Top: 
WXM. The data points show the count rate with 1 a uncertainty 
observed by detector XB, while the dotted line indicates the 
count rate observed by detector XA. We employ data from the 
full 2-25 keV band pass at a 6.6 s time resolution. The data is 
from the end of an observation, when the count rate decreases as 
the Earth covers an increasing part of the field of view. The flare 
is only visible in the XB light curve. From its position on the 
sky, 4U 1608-522 was in the field of view of XB and not of XA. 
Bottom: FREGATE light curve from the 6^0 keV band pass at 
1.2 s time resolution. 



lowed by a slow decay (Fig.|5]l. Since the event took place at the 
end of the observation, only one minute of the flare is recorded 
before the instruments are switched off. Nevertheless, the light 
curve suggests that the decay is much longer than for a normal 
type-I burst. The WXM detector image of the flare is consis- 
tent with the position of 4U 1608-522 if the spacecraft's attitude 
changed by 0.6° since the middle of the orbit. A variation of the 
attitude of this size is not uncommon at the end of an orbit. The 
FREGATE has no positional information, but 4U 1608-522 was 
in the field of view at that time. From both instruments spec- 
tral information is available in four energy bands. Due to the 
low number of observed photons per energy band, this informa- 
tion is of limited use for the WXM. For the FREGATE, the flare 
is observed in the 6 to 40 keV and 6 to 80 keV bands and not 
in higher energy bands, which is consistent with a type-I X-ray 
burst (such as a superburst). Comparing this flare to a type-I X- 
ray burst from 4U 1608-522 which was previously observed with 
the WXM, the flare's peak height is approximately 80% of that 
of the burst. Therefore, the peak flux of the flare is at most 80% 
of the Eddington limit. As the time of occurrence, the source po- 
sition and the peak flux are consistent with the superburst from 
4U 1608-522, we regard this as a likely observation of the start 
of this superburst. 

No precursor burst can be discerned. However, the presence 
of such a burst cannot be ruled out. If a potentially present pre- 
cursor is similar to the type-I bursts observed during outbursts, 
it would have a short e-folding decay time of a few seconds and 
a peak flux of around 60% of the Eddington limit. Before the 
superburst onset the statistical quality of the FREGATE data 
is such that we can exclude at a 3(7 confidence level a pre- 
cursor burst with a decay time of several seconds and a net 
peak flux of ~ 10% of the Eddington limit or higher However, 
the precursors observed from other hydrogen-accreting super- 



L. Keek et al.: First superburst from a classical low-mass X-ray binary transient 



7 



burster s all took place at the superburst onset dKuulkers et alj 
I2002bt IStrohmaver & M arkwardt"2002l: lin 't Zand et al.ll2003h . 
If the precursor and superburst light curves are superimposed, 
the presence of a precursor cannot be excluded from the data. 
Also, during the precursor burst the peak temperature might be 
lower than for the superburst, which would result in a lower ob- 
served flux in the FREGATE band. The WXM is sensitive down 
to lower energies, but the relative uncertainty in each data point 
is larger due to a smaller effective area, which may prevent us 
from detecting a precursor. 

After the superburst took place, the first normal type-I X-ray 
burst was observed with IBIS/ISGRI after 99.8 days. Since there 
are frequent data gaps, whose durations are long with respect to 
the duration of an X-ray burst, 99.8 days is an upper limit to the 
burst quenching time. 

In the 4.6 hours preceding the superburst (interval II, see Fig. 
lib) the persistent flux is at a level of 25.3cs^^ while in the 
7 days prior to that (interval I) it was 32.9 cs^'. The hardness 
ratio (Fig. 2};) indicates that the spectrum is harder in interval I. 
Approximately one day after the superburst peak, the count rate 
returns to the level of interval I. Only after 4.5 days does the 
count rate return to the level of interval II. This suggests that the 
superburst starts during a slight dip in the persistent flux. This 
makes it difficult to accurately separate the persistent emission 
from the burst emission during the decay. Therefore, we report 
results using both interval I and II. 

The decay of superbursts is typically well-fit by an exponen- 
tial. We fit the decay of the superburst with an exponential in all 
three channels separately as well as in the combined 1.5-12 keV 
energy band (see Table[3]and Fig.|4}l). In this fit we restrict our- 
selves to the first 5 hours of the decay, to avoid the tail where an 
assessment of the persistent emission is difficult. 

Not only is the count rate in interval II lower than in I. Also, 
the hardness ratio indicates that the spectrum is softer in interval 
II (see ). 

4.2. Spectral analysis, energetics and layer thickness 

We fit the persistent emission ASM spectra from both interval 
I and II simultaneously with an absorbed power law, coupling 
the A^H value for both spectra. The best fit values as well as the 
1.5-12 keV flux in each interval are provided in Table[3] 

We divide the superburst in four time intervals, as indicated 
in Fig.|4ll, from the part of the burst where most data points lie 
at > 3(7 above the persistent flux from interval II (when using 
interval I, we restrict ourselves to the first three burst intervals). 
These intervals contain in total 19 measurements. Note that some 
of these start at the same time, but are from different SSCs. A 
system atic error of 3%, de rived from the scatter in the Crab in- 
tensity dLevine et aljl996|). is in cluded in the uncertainties of the 
data points. ' Le vine et a" note that in many cases this is an 

underestimate of the systematic uncertainties. For each interval 
we extract a spectrum. We obtain the net-superburst spectrum by 
subtracting the persistent spectru m based on either inter val I or II 
from the observed spectrum (see lKuuIkers et al]|2002al for a dis- 
cussion on decoupling persistent and burst emission for normal 
type-I bursts). This net spectrum is fit with an absorbed black 
body. 

During the spectral analysis we fix the hydrogen absorption 
column density. From BeppoSAX broad-band observations dur- 
ing outburst we find a hydrogen absorption column density of 



Table 3. Superburst properties using persistent interval I and II 
with 1(7 uncertainties for the last digits between parentheses. 

Exponential decay times" 



Persistent: 


1 






11 






Band (keV) 


^exp 


(hr) 




'^exp 


(hr) 


xL/d.o.f. 


1.5-12 


4.11 


0.6 
0.5 


6.2/7 


4.8l 


-0.5 
0.4 


6.5/7 


1.5-3 


12.2 


+- 
-6.9 


1.6/6 


13.2 


+- 
-6.2 


1.6/6 


3-5 


7.81^:3 


1.4/7 


7.8l 


-2.3 
1.4 


1.5/7 


5-12 


o Q+0.5 
■^•■^-0.4 


2.7/7 


3.9^ 


-0.4 
0.3 


3.2/7 


Bolometric 


3.01 


U.4 

0.3 


4.3/7 


3.6l 


-0.3 
0.3 


5.7/7 



Persistent emission prior to flare (Nu coupled)* 


IntervaF 


I 11 






6.1(4)- 10^' cm"^ 




r 


1.75(3) 2.06(7) 




^powerlaw 


2.8(2) 3.3(4) 


photons keV 'cm ^s ' 


X^/d.o.f. 


1.7/1 




Flux'' 


1.3(3) 1.0(3) 


10^**ergs^' cm^^ 


Superburst spectral analysis (A'bb coupled)*"; persistent interval 1 


Interval"^ 


1 2 


3 


kT (keV) 


2.14(14) 1.90(12) 


1.60(8) 


A'bb 


1.4(3) -102 


xL'd-o.f. 




3.6/5 




3.18(6) 1.96(4) 


1.00(3) 


Superburst spectral analysis (N\,\y coupled)*"; persistent interval 11 


IntervaF 


1 2 


3 4 


kr (keV) 


2.13(12) 1.89(9) 


1.67(8) 1.38(7) 


A'bb 


1.7(3) -10^ 


xL'd-o.f. 




3.5/7 


Flux-^ 


3.81(2) 2.37(2) 


1.438(11) 0.68(6) 


Peak flux and burst fluence; persistent interval I 


Burst start: 


ASM HETE-2 


Maximum 


FluXpeak-' 


2.90(14) 3.5(3) 


4.8(6) 


Fluence ^ 


3.1(4) 3.8(6) 


5.2(1.0) 




1.5 1.7 


2.1 


En' 


1.4 1.6 


1.9 


Peak flux and burst fluence; persistent interval II 


Burst start: 


ASM HETE-2 


Maximum 


FluXpeak-' 


3.88(11) 4.5(2) 


6.0(5) 


Fluence ^ 


5.0(4) 5.9(5) 


7.6(1.1) 




3.0 3.3 


4.1 


En' 


1.5 1.6 


1.9 



A^H = 0.89 • lO^^cm- 
value of A^H = 0.61 • 



(see Sect. 13.21 ). However, here we use the 



"In the 1.5-3 keV band we find the decay to be consistent with a hori- 
zontal line, due to the large uncertainty of the data points in the tail of 
our fit interval. 

^Absorbed power law model (A'po^verlaw^ parameters: Nu the inter- 
stellar column density, F the photon index and N^ov/eilm/ the normaliza- 
tion defined as photons keV^'cm^^s^' at 1 keV. 
^^SeeFig.Elj. 

Unabsorbed flux in 1.5-12 keV energy range (lO^^ergs^' cm^^). 
*" Black body model parameters: temperature kr and normalization 
A'bb = -R^/rf^, with R the source radius in km and d the source distance 
in 10 kpc. 

■f^ Unabsorbed bolometric flux (10^** erss^' cm^^). 
* Unabsorbed bolometric fluence (10^ ergcm^^). 

Column depth (lO'^gcm^^). 
' Energy release per unit mass (lO'^ ergg^'). 



spectral analysis using similar ASM data and the same response 
matrix. 

In the spectral analysis, our best fit has a reduced X^ed value 
of 3.5. We increased the uncertainties in the data points such that 
we obtained a fit with ^ ^- Using this fit we determined the 

2 



10 cm which we obtained from the per- uncertainties in the fitted parameters using Ax = 1- Increasing 
sistent emission in interval II, since we perform the superburst the uncertainties to satisfy the "goodness of fit" condition Xred 
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was also required when we analyzed the Crab spectrum in Sect. 
12.21 In this procedure the best fit values of the model parameters 
deviate from the previous best fit values by at most 1 .37(7. We fit 
the four intervals simultaneously, keeping the value of the black 
body normalization the same for all spectra. Since the normaliza- 
tion parameter depends on the emission region size of the black 
body, this implies that the size is assumed not to change dur- 
ing the superburst decay. This is a valid approach, as we will see 
that the peak flux of the superburst is below the Eddington flux at 
which photospheric expansion is expected. Fits where we did not 
use this coupling did not indicate a significant change in radius 
from one interval to the next. We provide the best fit values of the 
model parameters in Table [3] Assuming isotropic radiation and 
using the distance of 3.2 ± 0.3 kpc as well as the best fit value 
for the black body normalization, we find an emission region size 
equivalent to a sphere with a radius ofR = 4.2 ±0.5 km, which is 
consistent wi th the b lack body radii both Mura kami et alj (Il980l) 
and Galloway et al.l (12006,) found for normal type I X-ray bursts. 

To investigate the influence of the choice of A'h on the fit 
results, we repeated the spectral analysis using A^h = 8.91 • 
10^' cm^^ (see Table[T]). The difference in best fit values of the 
model parameters to those in the previous analysis never exceeds 
1(7. The derived flux differs by at most 3%. 

Using the average count rate and the flux we obtained for 
each burst interval, we convert the count rate of each data point 
to flux. We fit an exponential function to the first 5 hours of 
the resulting curve and determine the peak flux at the start of 
this interval. We find a peak flux of 22 ± 3% (29 ± 3%) of the 
Eddington flux using the persistent level from interval I (II). 
Integrating this exponential gives us the burst fluence, which 
at a distance of 3.2 ± 0.3 kpc corresponds to a burst energy of 
(3.8±0.9)10'^ierg ((6.1 ± 1.2)10'*' erg). By extrapolating the 
exponential to the start of the flare observed by HETE-2, a peak 
flux of 27 ± 4% (34 ± 4%) of the Eddington flux and a burst 
fluence of (4.7 ±0.8)10'*' erg ((7.2 ±0.9) lO'*' erg) are obtained. 
Extrapolating to the start of the ASM data gap preceding the 
observed part of the burst, we find an upper limit to the peak 
flux and burst fluence: respectively 36 ± 6% (45 ±4%) of the 
Eddington flux and (6±2)10'*'erg ((9 ± 2)10^" erg) (see Table 

m. 

ICumming & MacbethI (|2004|) present multizone numerical 
models for the cooling of the neutron star surface layers dur- 
ing the decay of a superburst as well as analytic fits to these 
models. The models depend on the ignition column depth 
y = yi2l0'^gc m~^ and the energy release per gram E = 
/inlO'^ergg^'. ICumming et al.l (|2006 |) fit these models to the 
light curves of six superbursts. We perform the same fit to the 
light curve of the superburst from 4U 1608-522 using the per- 
sistent level from both interval I and II (Fig. |6] ). We find, re- 
spectively, yi2 = 1.5, £17 = 1.4 and yn — 3.0, £17 = 1.5. To 
estimate the values of yi2 and £17 in case the superburst started 
at either the start of the flare observed with HETE-2 or at the 
start of the data gap in the ASM light curve, we use the scal- 
ing laws found by Gumming et al. (2006). Thus we obtain the 
ranges yn = 1.5-2.1, £17 = 1.4- 1.9 and yn = 3.0-4.1, 
£17 = 1.5 — 1.9, using respectively the persistent level from in- 
terval I and II (see Table[3]l. 

5. Discussion 

The first detection of a superburst from a classical LMXB tran- 
sient poses significant challenges to superburst theory and un- 
derlines shortcomings thereof that were already realized earlier. 
This shows the need for additional ingredients in the models. 



Table 4. Gharacteristics of the 4U 1608 -522 supe r burst com- 
pared to 13 other superbursts (from i KuulkersI 12004 ' and 
in 't Za nd et al. 2004). Texp from fit to count rates converted to 
unabsorbed bolometric flux. 





Al J 1 /inO COT 


Other superbursts 


Duration (hr) 


5-22 


2-15 


Precursor 


NTn 


Always when data available 


texp(hl-)" 


~4.5 


0.7-6 


k7peak(keV) 


2.1±0.1 


1.8-3 


Lpe,k(103«ergs-1) 


0.4-0.7 


0.4-3.4 


ipers(-^Edd) 


0.11-0.15 


0.1-0.25* 


(Lpers)(iEdd) 


0.03 


0.11-0.43'' 


£b(1042erg) 


0.4-0.9 


0.3-1.4 


-r = £:b Apeak (hr) 


1.0-1.3 


1.1-6.9 


Y — ^ers/ ^peak 


0.2-0.5 


0.1-0.7 


^no bursts (d) 


<99.8 


>7.5 


Donor 


H/He 


He and H/He 



" Exponential decay time in 1.5-12 keV band for 4U 1608-522 and in 
similar energy bands for other superbursters. 



* Persistent luminosity at the time of the superburst, excluding 
the candidate ultracompact source 4U 0614+091 and the Z-source 
GX 1 7+2 which have a lumi nosity close to, respectively, 1% ^ Kuulkersl 
2005]; lin 't Zand et"aDl2007h and 100% (in 't Zand et al. 2004) of the 
Eddington luminosity. 

" Average persistent luminosity from all ASM observations, again ex- 
cluding 4U 0614+091 and GX 17+2. 




5 10 15 

Time (h) 



Figure 6. Fit of cooling models to superburst decay. The data 
points are ASM count rates converted to flux for each ASM mea- 
surement (see Sect. I4.2l i with 1(7 uncertainties. The persistent 
level from interval II is used (see Fig.|4]3). The solid and dotted 
line represent fits wi th, respectively, th e numerical and analytic 
cooling models from lGumming & Macbeth (.2004) . 



While the superburst characteristics of 4U 1608-522 are similar 
to those from other superbursters (see TablelDi, the accretion rate 
on a time scale of ten years is a factor of 3 smaller than observed 
for the slowest accretor among all other superbursting sources. 

In the following subsections we will see that there are two 
challenges: the destruction of carbon by the ever-present hydro- 
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occurred. This implies that the carbon fuel has been produced 
during multiple outbursts. 



5.2. Producing and preserving carbon 



50500 51000 51500 52000 52500 53000 53500 54000 
Time [MJD] 

Figure?. Accreted column depth yi2 = ^/(lO'^gcm^^) over 
time. Assuming a neutron star of 1 .4M0 mass and 10 km radius 
at a distance of 3.2 kpc, we show the cumulative column depth 
that is accreted since the start of RXTE ASM observations as in- 
ferred from the X-ray flux at 14 day time resolution (c.f. Fig.[T]). 
The lower curve takes into account only matter that is accreted 
when 4U 1608-522 was in the high state when the ASM count 
rate exceeded 20 c s^^ while the upper curve includes the low 
state as well. The dotted line indicates the moment the super- 
burst was observed. 



gen/helium flashes is apparently not as important as theory pre- 
dicts and the predicted crust temperatures fall short by a factor 
of 2 of the value required for the ignition of unstable carbon 
burning. Additionally, we make a few remarks on the mHz QPO 
observed in 4U 1608-522. 



5.1. Recurrence time 

From fits of cooling models to the superburst light curve we 
find that the superbursting layer had at ignition a column depth 
of yn — 1.5—4.1 and an energy release per unit mass of 
£17 = 1.5 — 1.9, which implies a total energy release of (2.6 — 
9.8)10'^2gj.g From the X -ray flux we derived an energy release 
of (0.4 - 0.9)10'*2ej.g, which is only 5 - 45% of the total in- 
ferred from the cooling models (after correction for the gravita- 
tional redshift with a typical value of 1 +z = 1.3). The rest of the 
energy is either released as neutrinos or is conducted to greater 
depths and radiated away on a longer timescale. Fig.|7]shows the 
column depth that is observed to be accumulated since the start 
of the ASM observations. At the moment of the superburst a col- 
umn of = 0.55 (^12 — 0.27 taking into account only the stable 
burning in the high state) had been accreted over a period of 9.3 
years. Note that this is derived under the assumption that the X- 
ray flux is a good measure for the mass accretion rate, which 
may not be the case (see discussion in Sect. |3.2| i. In fact, we as- 
sume that stable burning occurs when the ASM flux is in excess 
of 20 c s^', which introduces an uncertainty of several tens of 
percents in the column depth. With an average accretion rate of 
M = ( 1 .8 ± 0.2) lO^^Medd the required column depth is accreted 
over a period of 26 to 72 years. If carbon is only produced from 
matter accreted during stable burning in the high state, the recur- 
rence time is twice as long. In both cases it is much longer than 
the 57.6 day duration of the 2005 outburst before the superburst 



ISchatz et al.l (l2003h found from nucleosynthesis calculations 
that carbon is produced most efficiently through the stable nu- 
clear burning of hydrogen and helium. Unstable burning dur- 
ing type-I X -ray bursts results in significantly less carbon (< 
1% by mass; Woosley et akl 2004 ). but produces heavier ashes. 
ICumming & BildstenI (12001 ) proposed that these heavy ashes 
are important because they set the opacity and the refore tem- 
perature of the superbursting layer, althoughlBrownl (2000) later 
showed that this effect is much smaller than initially estimated 
once the temperature profile of the entire crust was taken into 
account. Therefore, it is interesting to ask whether both burn- 
ing regimes are necessary for producing the fuel for superbursts. 
Certainly, all known superbursters exhibit both type-I X-ray 
bursts and a high a value, which are indications for both un- 
stable and stable burning, respectively. 

The stable burning of helium takes place at accretion rates 
M in excess of 10% of Eddington. Therefore, we expect that in 
4U 1608-522 carbon is created during the outbursts. This carbon 
could be destroyed by type-I X-ray bursts, which are observed 
during all flux states. For carbon to survive the bursts, it has to 
reside in a layer which is depleted of hydrogen. IWoosley et al.l 
(l2004h find this to be the case at the bottom of a bursting 
layer The different composition of the ashes of unstable hy- 
drogen/h elium burning as opposed to sta ble hydrogen/helium 
burning ( Woosley et"al1l2004tlSchatz et al.|[ 2003) will likely mix 
rapidly. If a period of stable burning was followed by unstable 
burning, the heavier ashes of the latter would rapidly undergo a 
Rayleigh-Taylor instability and mix with the lighter ashes below. 

An important question is whether the carbon i s able to burn 
stably during the low state. ' Gumming & BildstenI (1200 Ih argued 
that carbon would burn stably at low accretion rates, and sug- 
gested that this may explain why superburst sources have accre- 
tion rates of « 0.1 Eddington and up. This question is crucial to 
understand for 4U 1608-522, since the accretion rate is low dur- 
ing quiescence and the low persistent state. Further modelling 
of the production and destruction of carbon during the cycle of 
transient outbursts is needed. 

A complication here is that at the low accretion rates that oc- 
cur during quiescence and in the low persistent state in 4U 1608- 
522, the accretion may be slow enough that the carbon and 
heavy elements in the ashes of normal type-I bursts gravitation- 
ally separate, with the carbon floating upwards relative to the 
heavy elements. To estimate the relative drift velocity of a car- 
bon nucleus in the heavy element ocean, we follow the work 
of Bildsten & Hall (2001), who estimated the rate of sedimenta- 
tion of ^^Ne in the interiors of CO white dwarfs. They start with 
the self-diffusion coefficient for a one component plasma, D w 
ScOpG^T^'^/^ dHansen et al. I ll975h . where (Op is the ion plasma 
frequency, a is the inter-ion separation, and F is the standard 
parameter that measures the strength of the Coulomb interac- 
tions between ions. A single carbon nucleus in a background 
of ^^Fe feels an upward force of 12mp^/13, giving a drift ve- 
locity V — 54mpg/169eri/^(4;rp)i/2 (compare Bild sten & Halll 
120011 eq. [3]). Writing the time to cross a pressure scale height 
at this drift velocity as = H/v = P / pgv, and assuming that the 
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pressure is provided by degenerate relativistic electrons (a good 
approximation at these depths), we find 
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where we assume that the heavy element is ^^Fe (the scaling is 
ts °= Z^^/'^^/A). This timescale is equal to the accretion timescale 
facer for m K, 0.05 medd- The sedimentation timescale ts is com- 
parable to the timescale on which we infer the carbon fuel is 
accumulated, so that the physics of the relative separation of the 
carbon and heavy elements should be considered further. 

5.3. The crust temperature and carbon ignition 

Although the accretion rate of 4U 1608-522 reaches M ^ 
0.23MEdd during the peak of the outburst, the high accretion 
rate only lasts a short (Ss 100 d) time, so that the crust does 
not reach a thermal steady-state. The temperature of the deep 
crust can be inferred by observations of the quiescent luminos- 
ity (Rutledge et al. 1999). Hence this source provides a test of 
the strength and location of heat sources in the crust. The ques- 
tion is whether the crust can reach temperatures sufficiently hot 
to initiate thermally unstable '^C fusion in the short time avail- 
able during the outburst. To check this, we constructed a time- 
dependent model based on recent calculati ons of the heatin g 
from electron capture reactions in the crust jGupta et al.ll2007l) . 
We set the mass number to A = 40, which yields a large heat 
release (0.40 MeV per accreted nucle on) in the outer crust. Fo r 
the inner crust, we use the model of Hae nsel & ZdunikI (Il990h . 
We integrated the thermal diffusion equation using a standard 
m ethod-of-lines form alism. The microphysics is identical to that 
in lGupta et alj (l2007l) . The total heat released in the crust during 
outburst is LSMeVu"'. 

As an initial condition, we computed the temperature for 
steady accretion at the time-averaged mass accretion rate (M) ~ 
6.5 • lO'^ gs^'. We then ran through a series of 60 out- 
burst/quiescent cycles, to ensure that the crust had reached a 
limit cycle. During outburst we set M — 1.3 ■ lO'^ gs^', with 
a duration of 80 d; the quiescent interval was set to 4.16 yr. 
During quiescence we used a relation between the temperature 
at a column 3' = 10^ gcm^^ an d the surface lumin osity appro- 
priate for an accreted envelope ( Bro wn et al.l 120021) . The quies- 
cent luminosity reached a minimum ofLg = 6.1- 10^3 erg s ' , in 
reasonable agreement with the value derived by iRutledge et alj 
(I2OOO1) . Figure [8] shows the change in temperature, as a func- 
tion of column, during the outburst. At the start of the outburst 
(solid line), we set the temperature at a column y = 10^ gcm^^ 
to 5 • 10** K. This is a reasonable upper limit on the temperature 
set by steady-state H/He burning. After 80 d, the temperature in 
the outer crust has risen substantially (dotted line), but is still 
much cooler than needed to explain the fitted ignition column 
and temperature (shaded box). 

We conclude that heating in the crust from electron captures, 
neutron emissions, and pycnonuclear reactions, is insufficient to 
raise the crust temperature during the outburst enough for to 
ignite at the inferred column. 

In the other superburst sources, which are either persis- 
tent accretors, or transients with long duration outbursts (e.g. 
KS 1731-260), the crust reaches a thermal steady state. Even 
in these sources, theoretical models have difficulty explain- 
ing the inferred su perburst ignition coliimns of ~ lO'^ g cm~^ 
dBrownl 12 004; Co oper & NaravanI 120051: ICumming et aL[|2006t 
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Figures. Temperature evolution in the crust of 4U 1608- 
522 during the outburst. At the beginning of the outburst (solid 
line) we set the temperature at a column y = lO^gcm^^ to 
5 • lO^K. After 80 days, the temperature in the outer crust has 
risen (dotted line), but is still well below the temperature for 
which a one-zone stability analysis would predict a superburst 
with a depth inferred from fits to the cooling (shaded box). 



core that affect the thermal profile such as neutrino emissivity 
and crust conductivity must all be chosen to maximize the crust 
temperature in order to achieve the observed ignition depths. In 
fact. lCumming et alj (l2006l) showed that it was impossible to re- 
produce the observations if neutrino cooling due to Cooper pair 
formation of superfluid neutrons in the crust (e.g. lYakovlev et al] 
1999) was included, because this efficient neutrino emission 
mechanism led to an upper limit on the crust temperature at 
the ignition depth of « 4 x 10^ K, less than the temperature 
of » 6 X 10^ K required for carbon ignition. Recent work by 
iLeinson & PerezI (2006) has shown that the neutrino emission 
due to this mechanism has been overestimated in previous cal- 
culations, and that this mechanism is in fact not the dominant 
source of neutrino emissivity. 

As we have shown, the short duration of the transient out- 
bursts in 4U 1608-522 means that the crust is heated to much 
lower temperatures than those corresponding to thermal equilib- 
rium. Therefore the observation of a superburst from 4U 1608- 
522 is deeply puzzling and presents a significant challenge to 
current superburst models. 



5.4. Bursting behavior and mHz QPOs 

The normal type-1 X-ray bursting behavior looks similar to 
what is seen in persistent bursting sources with strongly variable 
persistent fluxes, most notably 4U 1820-30 3 ( Chou & Grindlavi 
2001), GX 3-H l (Iden Hartoe et all j2003l) . 4U 1705-44 and 
KS 1731-260 dCornehsse et al.1 l2003h . In those cases X-ray 
bursts become far less frequent or even absent during high 
flux states. Also in the case of 4U 1608-522 the bursting be- 
havior is seen to diffe r between the low and high flux state. 
[Murakami et al.l (Il980l) report that in the high flux state bursts 
have a shorter duration and a higher a value than the bursts in the 
low flux state, while the burst rate is not significantly different. 
The bursts observed with the RXTE PCA and the BeppoSAX 
WFCs exhibit this behavior as well. This would be consistent 
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with hydrogen and heUum burning becoming predominantl y sta- 
ble at higher mass accretion rates ("van Paradiis et alj|198 8'). 

fRevnivtsev et al. (2001) examined EXOSAT and RXTE 
PCA data of 4U 1608-522 and discovered low frequency quasi- 
periodic oscillations (mHz QPOs). They observed no type-I X- 
ray bursts when the persistent emission was higher than when 
the QPOs were detected, which supported the hypothesis that 
mHz QPOs occur at the transition of unstable to stable burning, 
such that at higher accretio n rates the bursting behavior ceases 
(see also iHeger et al.ll2007l) . Interestingly, the PCA, WFCs and 
IBIS/ISGRI have observed a total of four X-ray bursts when the 
persistent flux was up to a factor of two higher than when mHz 
QPOs were seen (see Fig. [3]). Although this disagr ees with the 
prediction of the marginally stable burning model of lHeger et al.l 
( 20071). the hypothesis could still be valid if, as discussed by 
iHeger et aTl (|2007), the conditions for nuclear burning vary sub- 
stantially across the neutron star surface, such that on one part of 
the surface unstable burning gives rise to type-I bursts while on 
another part mHz QPOs occur. The scenario that ignition con- 
ditions might vary across the surface of the star for rap idly ro- 
tating neutron stars has been discussed by Cooper & N aravanI 
( l2007h (see also lSpitkovskv et al.ll2002l) . although the predicted 
range for the global accretion rate where both stable and unsta- 
ble burning take place at different latitudes is a factor of 3 higher 
than where mHz QPOs are observed from 4U 1608-522. 

6. Summary and conclusions 

For the first time a superburst has been observed from a tran- 
siently accreting neutron star in a low-mass X-ray binary where 
the duration of the accretion outbursts is short with respect to 
the expected superburst recurrence time. 4U 1608-522 exhibited 
a superburst 55 days after the onset of an accretion outburst. We 
analyzed the superburst as well as the long term accretion and 
bursting behavior of 4U 1608-522. We found the properties of 
the superburst to be comparable to those of the long superbursts 
observed from 4U 1254-690 and KS 1731-260. Fits of cooling 
models to the superburst light curve indicate that the superburst 
ignited when a column of y = (1.5 — 4.1)10^^gcm^^ had been 
accreted. From the average X-ray flux that the RXTE ASM ob- 
served from 4U 1608-522 over the past 1 1 years, we derive that 
this column was accreted over a period of 26 to 72 years. 

The difficulties that superburst theory has in explaining the 
ignition column depths derived for other superbursters, are un- 
derlined by the observation of the superburst in 4U 1608-522. 
The transient nature of the accretion poses challenges for super- 
burst models. The average accretion rate over the past 1 1 years 
was much lower than for the oth er superbtirsters , with the possi- 
ble exception of 4U 0614H-091 ( Kuulkerslllool) . At the derived 
column depth, current superburst theory predicts that during the 
timespan of the outbursts the crust temperature does not rise to 
the value required for superburst ignition. Furthermore, the car- 
bon that fuels the superburst is expected to be produced only 
during the outbursts, while the stable nuclear burning of ^^C as 
well as the frequent type-I X-ray bursts lower the carbon abun- 
dance. The detection of the superburst implies that a significant 
amount of carbon survives the long periods outside the outbursts. 
The balance between the creation and destruction processes is 
influenced by mixing and sedimentation in the neutron star enve- 
lope. The effect of these requires further study before predictions 
about the produced amount of superburst fuel can be made. 
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